Local atom number fluctuations in quantum gases at flnite temperature 
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We investigate the number fluctuations in small cells of quantum gases pointing out important 
deviations from the thermodynamic limit fixed by the isothermal compressibility. Both quantum 
and thermal fluctuations in weakly as well as highly compressible fluids are considered. For the 2D 
superfluid Bose gas we find a significant quenching of fluctuations with respect to the thermodynamic 
limit, in agreement with recent experimental findings. An enhancement of the thermal fluctuations 
is instead predicted for the 2D dipolar superfluid Bose gas, which becomes dramatic when the size of 
the sample cell is of the order of the wavelength of the rotonic excitation induced by the interaction. 

PACS numbers: 67.85.-d,03.75.Hh, 42.50.Lc, 05.30.Fk, 05.40.-a 



I. INTRODUCTION 

The experimental possibility of detecting few and even 
single atoms confined in magnetic or optical traps with 
high precision is opening new perspectives in the study 
of correlations i2j] and fluctuations in atomic and 
molecular gases at the mesoscopic and microscopic scale. 

In the present work we investigate the fluctuations of 
the particle number in small sample cells of quantum 
gases, exploring the deviations from the predictions of 
thermodynamic theory which relates the fluctuations to 
the isothermal compressibility. 

We consider two important examples of continuous sys- 
tems. The first one is the ideal Fermi g benchmark 
of statistical mechanics, which has been recently stud- 
ied experimentally in atomic gases to point out the anti- 
bunching effect produced by the Pauli exclusion principle 
d, Q . The second system is the dilute quasi- two dimen- 
sional (2D) bosonic gas in the superfluid regime. Two 
dimensional gases are particularly suited to measure the 
number fluctuations since one avoids the column density 
integration. Measurements of this kind have been re- 
cently carried out showing important quenching of the 
fluctuations in the degenerate regimeQ, whose origin is 
clearly explained by our analysis. On the other hand 
quasi-2D Bose gases with dipolar interaction present an 
enhancement of number fluctuations. This is due to the 
roton-like spectrum which has important conse- 

quences on the behaviour of the static structure factor, 
thus on the number fluctuations, especially at finite tem- 
perature. 

When the size of the sample cell is large enough one can 
use the thermodynamic theory of fluctuations (see, e.g., 
(m ) . The expression for particle number fluctuations at 
temperature T reads 



NkBTriKT, 



(1) 



which relates the fluctuations to the isothermal compress- 
ibility kt = l/n^dn/dpL = 1/nmc^, with n the average 
density, /i the chemical potential and c the speed of sound 
[l^ . When T is large the compressibility of the gas ap- 
proaches the classical value tikt = i/ksT and one recov- 



ers the shot noise regime. An extension of the thermody- 
namic relation for the density fluctuations of non-uniform 
systems has been recently proposed as a tool to measure 
temperatures If the radius of the cell is not large 

enough, important deviations with respect to Eq. ([T]) 
start occurring. In our investigation we will mainly focus 
on the low temperature regime where quantum effects 
become particularly important. 

The paper is organized as follows. In the next section 
we recall how to calculate the particle number fluctua- 
tions from the correlation function and the static struc- 
ture factor, respectively. In section III we discuss the 
temperature dependence of fluctuations in a flnite volume 
of an ideal Fermi gas in three dimensions. In the same 
section we also study the dependence of quantum fluctua- 
tions (T — 0) on the shape of the probe volume. Then in 
section IV we analyse the density fluctuations for the case 
of an interacting quasi-2D Bose gas at temperatures far 
below the critical temperature. We condsider short-range 
interactions as well as long-range dipole-dipole interac- 
tions and discuss differences in the static structure factor 
and in the flnite volume fluctuations. In the same sec- 
tion we provide an analytical expression - whose deriva- 
tion for completeness is given in the Appendix - for the 
zero-temperature fluctuations in 2D and for a contact 
potential. 



II. FLUCTUATIONS IN A FINITE VOLUME 



The particle number fluctuations {SN'^) in a flnite 
cell of volume V are calculated by double integration of 
the density-density correlation function nv(\ri — 
(ri(ri)n(r2)) — over V and can be written as[l^. n 



(SN^) = n v{r)h{r)dr = n S{k)H{k) 



dk 



i2n) 



D ■ 



(2) 



where D is the dimensionality of the system, h{r) = 
Iv Iv ^(''i — r2 — r)(iri(ir2 is a geometrical factor which 
depends on the shape of the probe cell, H{k) — 
/y /y ^*''^'^^~'"^'''^''i'^'*2 is its Fourier transform and 



2 



S'(k) = J i^(r)e~''""dr is the static structure factor. Geo- 
metrically h{r) corresponds to the volume of the over- 
lapping region between a cell of volume V and the 
same cell shifted by the vector r. For large volumes 
(thermodynamic limit) i?(k) — ^ V6{'k) and one finds 
{6N'^) — NS{k — 0) yielding, at finite temperature, the 
thermodynamic expression Eq.([T]). Deviations from the 
thermodynamic limit are directly related to the finite- 
k behaviour of the static structure factor. In order to 
discuss these deviations it is useful to distinguish be- 
tween weakly compressible and highly compressible flu- 
ids. Weakly compressible fluids are characterized by the 
condition mc^ > Edeg where Edeg ~ fi?n?/^ /m is the de- 
generacy energy and include systems like the ideal Fermi 
gas, the unitary Fermi gas and superfluid Helium. For 
such systems the low temperature condition fc^T ^ mc^ 
is easily reachable, yielding antibunching with respect 
to the shot noise limit, i.e., SN"^ < N. Highly com- 
pressible fluids are instead characterized by the condi- 
tion mc^ ^ Edeg- For such systems, which include di- 
lute superfluid Bose gases, the low temperature condition 
ksT <C mc^ is hardly reachable experimentally and the 
fluctuations naturally exhibit bunching effects with re- 
spect to the shot-noise limit. In the following we will 
consider two important examples: the ideal Fermi gas 
and the 2D Bose gas. 



III. IDEAL FERMI GAS 

As a paradigmatic example of weakly compressible flu- 
ids we consider the ideal Fermi gas, whose flnite tempera- 
ture density-density correlation function in three dimen- 
sions reads 



v{v)^5{v)- 



ksT kl 
nep 27r^r 



oo u sm 
du 



in {u\J 



kpr 



(3) 



where kp = (67r^n)^/'^ and ep — h'^kjp/2m are the Fermi 
momentum and the Fermi energy, respectively. The 
chemical potential of the ideal Fermi gas /x(T) is obtained 
self-consistently from 



keT 



3/2 



du 



3u^ 



= 1 



(4) 



the larger the volume (and hence N) and, for finite vol- 
umes, they approach each other at large temperature 
when the deBroglie wavelength becomes smaller than the 
sample size. In particular the fluctuations in a flnite vol- 
ume are always larger then the ones predicted by the 
thermodynamic expression Eq. ([l} and remain finite at 
T = where they have a purely quantum nature. In this 
limit, to the leading order in N, they can be written 
as (B^nl (SN^) « 3(3/(327r4))i/37v2/3ijj(io.45iVi/3). 
Such an expression suggests that quantum fluctuations 
have a surface nature being dominated by the depen- 
dence iV^/^ c>c R^. Let us also note that the approxima- 
tion (SN^) « nNkBHT=o = 3NkBT/{2ep) obtained by 
neglecting the T-dependence in the compressibility (see 
dashed-dotted line in Fig. [T]), is bad for all T: at low 
temperatures quantum fluctuations are not taken into ac- 
count and at large temperature the relation kt ~ ht=o 
is not valid. 



A. Surface shape dependence of quantum 
fluctuations 

In order to get a better insight into the geometrical de- 
pendence of the quantum fluctuations [1^ we calculate 
((5iV^) at T = for a cylinder with radius R and length 
L, keeping the volume (and hence N) flxed. The par- 
ticle number fluctuations are given by Eq. ([2]), with the 
geometrical factor 



h(p,z) = {2L-z) 



2i?2 arccos f — ) - - \/^R^~7^ 
\2R/ 2 



(6) 



where z G [0, 2L] and p G [0, 2R]. The results are shown 
in the inset of Fig. [1] as a function of the aspect ra- 
tio L/R and for N = 100. We flnd that the particle 
number fluctuations are minimal for L R where the 
surface is minimal and they are even smaller for a sphere 
with the same volume (horizontal dashed-dotted line). 
Interestingly, in the inset of Fig. [T] we also show that 
the result follows closely the dependence of the surface 
Scyi = 21/2/3/^1/3 [2(i/i?)i/3 + (i?/L)2/3] of the cylin- 
der on the ratio L/R, for a flxed volume V . Deviations 
from the latter expression is due to the logarithmic term, 
which cannot be taken into account by the geometrical 
analysis. 



Initially, we consider the particle number fluctuations 
in a spherical cell of radius i?, for which the geometrical 
factor is 



AttR-^ 



- 7:rR^ 



TTr 
12 



(5) 



with r e [0,2i?]. In Fig. [T] we show (SN'^) as a func- 
tion of temperature for N = A/SnR^n = 100 and 2000 
calculated from Eq. (solid lines) and from Eq. ([T]) 
(dashed line). As expected the two results are the closer 



IV. 2D BOSE GAS WITH DIPOLAR 
INTERACTION 

As important examples of highly compressible fluids 
we analyse dilute bosonic gases in the superfluid regime. 
Moreover we consider the system to be in a quasi-2D ge- 
ometry where the motion along z is frozen. This is usu- 
ally realised with a strong harmonic conflnement in the 
z direction. We consider a quasi-2D gas, where the three 
dimensional s-wave scattering length a is much smaller 
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dispersion relation e{k) of the elementary excitations via 



FIG. 1: Number fluctuation (SN^) for a free Fermi 
function of T/Tf calculated integrating the density-density 
correlations (sohd line) for A'^ = 100 and A'^ = 2000 and using 
the thermodynamic expression Eq. ^ (dashed line). Inset: 
{SN^) at T = for a cylinder (solid line) and for a sphere 
(dashed-dotted line) for A'^ — 100 as a function of the aspect 
ratio L/R. The fluctuations closely follow the cylinder surface 
dependence (dashed line) for a flxed volume (see text). 



than the harmonic oscillator length Iz in the z direction. 
In this case it is possible to define an effective 2D short- 
range coupling constant g = h^\/%T:a/ [mlz). We will 
also include a long range dipole-dipole interaction, char- 
acterised by the parameter ga — 2y/STT(P / (3^z), where d is 
the dipole moment of the atom (or the molecule) which 
we take oriented perpendicular to the two dimensional 
plane of the gas. 

Let us introduce the healing length ^ = h/{mc) and 
the chemical potential /i — mc^ — h^/{m£,'^) = ^(1 + 
f3)n, where /? — ga/g — 2m(P / {?>It? a) is the ratio be- 
tween of the dipole-dipole and the contact interaction 
strengths f2Q|. Within Bogoliubov theory the spectrum 
of elementary excitations is [l^l 



e(fc) 



1 + /3 20F 



F 



klz_ 
V2 



(7) 



with F{x) = xerfc (x) . In the absence of dipolar in- 
teraction {f3 — 0) this spectrum has the well-known Bo- 
goliubov form. For finite /3 the dipolar term can signif- 
icantly change the shape of the spectrum. For negative 
s-wave scattering length a the spectrum can exhibit a 
roton-like minimum. Note, that in this case the gas must 
be stabilised against phonon instability, which demands 
1/31 >1. 



5(fc,T) 



2l2 



2me(fc) 



■ coth 



(8) 



The latter relation, Eq. (|8]), follows from the fact that 
the imaginary part of the response function in a weakly 
interacting Bose gas is T-independent and takes the form 
of a (5-function 19f]. 

The static structure factor depends significantly on the 
ratio UbT / {mc^). In Fig. [5]we show S'(fc, T) as a function 
of fc^ for different temperature with and without dipolar 
interaction. The peculiar roton-like shape of e(fc) shows 
up in a maximum in S{k). This effect is significantly 
amplified at finite temperature as can be seen in Fig. [21 
Actually, even if the spectrum does not exhibt a roton 
minimum (see the inset in the figure), a clear peak is 
visible in the static structure factor. Physically, this is 
due to thermal occupation of the roton-like states at finite 
k. The enhancement can be easily understood from Eq. 
dl]) since, if e{k) <C fcsT, the static form factor reduces 
to 



S{k,T) ~ S{k,Q) 



2kBT _ tfk^kBT 
e(fc) me(k)'^ ' 



(9) 



and in particular at low momenta we have S{k^ T) ~ (1 + 
AkC)kBT/{mc^) for ga ^ 0, with A = 3/3;^/[v^(l + /3)] 
and S{k,T) = (1 - {k£, / 2f)k bT / {mc^) for ga = 0. It 
is worth noticing that the large thermal enhancement of 
the the structure factor near the roton minimimum is pe- 
culiar of dilute gases. In fact in a weakly compressible 
system, like superfluid helium the thermal energy fc^T is 
always smaller than the roton energy. Moreover Bragg 
spectroscopy, a technique typically used to extract S{k) 
in quantum gases, can only access the T = static struc- 
ture factor, being sensitive to the imaginary part of the 
response function, rather than to the dynamic structure 
factor 0. We will show that the study of the particle 
number fluctuations can instead reveal the effect of the 
sizable temperature enhancement in S{k) produced by 
the rotonic excitation. 



B. Fluctuations without dipolar interactions 

First, we study the temperature dependence of the 
fluctuations in the non-dipolar case {ga — 0). In order 
to calculate the particle number fluctuations in a disk of 
radius R we use Eq. ([2]) with the geometrical factor 



H2D{k) = 



An^R^jfiRk) 



(10) 



Static structure factor 



At low temperatures T <^T\, where the gas is super- 
fluid, the static structure factor S{k,T) is related to the 



where Ji denotes the Bessel function. At T = the 
fluctuations can be calculated analytically (see appendix 
A) 



mg ^ 



(11) 
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FIG. 2: Static structure factor S{k, T) for a quasi-2D Bose 
gas with (continuous lines) and without (dashed hues) dipolar 
interaction for kBT/mc? — 0, 1, 2. Inset: the dispersion 
relation Eq. 0. For the dipolar case we use values for ^^Cr 
atoms, namely a — — 1.5ao, Iz = 0.09/^m which corresponds 
to /3 = -1.083, n/{V2^lz) = 10^'' cm"^ and ^ = 2.5^m. 



with C w 7.4. The result for SN'^ as a function of tem- 
perature and for N — 12, 100 is shown in Fig. [3]where we 
also report the result ([1]), which holds in the thermody- 
namic limit. We used the parameters of the experiment 
reported in Ref. In the inset we plot the fluctua- 

tions versus the interaction strength and compare with 
values measured in the superfluid regime of experiment 
Q • The agreement is pretty good. Further we show that 
the fluctuations are closer to the T = result Eq. pT|) 
(dashed-dotted line) than to the thermodynamic limit. 
At finite T we calculate the fluctuations using Eq. (|5]) 
with the dispersion relation Eq. ([7]) ioi (3 — 0. 

Fig. [3] shows that we have a significant deviation from 
the thermodynamic relation Eq. ([T]) at all temperatures. 
We can explain the behaviour of the fluctuation in an 
easy way within a simple approximation which describes 
how the fluctuations are related to the size of the sample 
cell, indeed if we consider a sphere of radius R the rele- 
vant momenta contributing to the integral ([2]) are of the 
order and the fluctuations can be approximated by 
the expression SN'^ ~ NS{a/R, T), where a is a constant 
of order unity. There are three natural length scales in 
the problem: the healing length ^, the phonon thermal 
wavelength Ip = hc/ksT and the de Broglie wavelength 
It — h/\/'mkBT. 

If i? S> ^ the static structure factor reduces to 



S ki 



kBT\ 



( rru? ai^ 

coth 

2i? yiksT R 



^cothf^^ 
2i? \2R 

(12) 

and, as expected, only the phonon part of the spectrum 
plays a role. If the radius is also larger than the phonon 
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FIG. 3: {SN^)/N for a quasi-2D Bose gas in the superfluid 
regime as a function of ksT / {m(?) for TV — 12, 100 and g — 
/4m, calculated from Eq. ((2|) (solid lines) and from Eq. ([T]) 
(dashed lines). Inset: {5N'^) as a function of the inverse of the 
dimensionless coupling constant g = mg/f? for A'^ = 12, T — 
2QnK calculated from Eq. ^ (solid lines) and from Eq.([T]) 
(dashed lines). The (red) small boxes are the experimental 
value extracted from [tI at T 20nK. The T = resuh Eq. 
(Illf) is also reported for comparison (dashed-dotted line). 



thermal wavelength [R ^ Ip) we get, as expected jll|, 
SN"^ = NksT /m(? , i.e., we recover the thermodynamic 
result and, depending on the value of T /mc^, the fluctua- 
tions can be sub- or super-poissonian. If, instead, R <^ Ip 
we have SN^ = aNS,/2R cx N^^^, i.e., aside from the log- 
term, we recover the T = quantum result. Notice that 
the latter regime requires the condition ksT <^ mc? and 
hence is difficult to access in dilute Bose gases. 

In the opposite i? <C ^ regime, we expect to probe 
the particle-like part of the spectrum. Indeed the static 
structure factor reads 



S 



coth 



mc 



:coth 



' (13) 

Again we can distinguish two cases: for R ^ It, 
which is possible only for ksT ^ mc^, we get JTV^ = 
NkBT/mc^{R/£)'^, i.e., the fluctuations are reduced with 
respect to the thermodynamic limit, as clearly shown in 
Fig. [3] for kBT/mc^ > 1. If instead R < h, i.e. if R is 
the smallest length scale of the problem, we recover the 
shot-noise result SN"^ = N which, by the way, coincides 
with the T = ideal Bose gas result. 
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FIG. 4: {5N^)/N as a function of R/(^ at zero temperature 
(dashed lines) and ksT / (mc?) = 2 (solid lines) for a dipolar 
gas of ^^Cr atoms (black lines) and for a non-dipolar gas (/3 = 
0) (grey lines) with the same chemical potential. Parameters 
as in Fig. [2] and in particular ^ = 2.5^m. 



C. Fluctuations with dipolar interactions 

In order to calculate the particle number fluctuations 
in the presence of dipolar interactions {gd 7^ 0) in a disk 
of radius R we use again Eq. ^ with the geometrical 
factor [ini In the previous subsection we have seen that 
a discussion of the proper length scales can give a good 
estimation for the behaviour of the flutuations. For a 
dipolar gas there is another important length scale in the 
problem: 1/kmax, which corresponds to the maximum of 
S{k). Notice that even a small deviation of the excita- 
tion spectrum from the usual Bogoliubov form, without 
exhibiting a minimum (see inset in Fig. [2]) leads to the 
maximum in S{k), which is strongly amplified at finite 
T as shown by Eq. Hence, for sample cells of size 

R ^ we expect, at finite T, a significant am- 

plification of the particle number fluctuations. In Fig. 
lU we report SN^ calculated for a dipolar gas of ^^Cr 
atoms under experimentally feasible conditions [l^. No- 
tice that the amplification of the fluctuations survive also 
for larger values of the cell size and the measurement of 
the particle number fluctuations is expected to provide 
a very sensitive tool to reveal the temperature amplifi- 
cation of the static structure factor in the rotonic-like 
region. 



V. CONCLUSION 

We have shown that the investigation of the atom num- 
ber fluctuations at the local scale provides a new insight 
on the microscopic structure of the correlations present 
in quantum gases where quantum and thermal effects 
combine in a non-trivial way. In particular, we explain 
the quenching of atom number fluctuations measured in 
a quasi-2D superfluid Bose gas Q- In a quasi-2D dipolar 
gas we obtain a crucial thermal enhancement of fluctua- 



tions. This opens new perspectives for future experimen- 
tal investigations. The present analysis can be naturally 
extended to investigate also spin-fluctuations in quantum 
systems [Hill]. 

A more quantitative comparison with experiments, 
where the cell V has not sharp boundaries, would re- 
quire the inclusion of the proper smooth weight function 
in the convolution integral Eq. ^ . 
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Appendix A: Zero temperature fluctuations of the 
2D Bose- gas 

In the present appendix we calculate analytically the 
particle number fluctuations at T = from Eq. ^ for a 
disk with radius R. Since the spectrum is given by Eq. 
([7]) for /3 = 0, and the geometrical factor by Eq. (fTO|) . the 
integral giving the number fluctuations can be written as 




2 4 6 8 



R/4 

FIG. 5: Comparison of (^iV^) as a function of R/^ for g — l/A 
calculated numerically from Eq. (|A1|) and Eq. HA7|) . Inset: 
Zoom onto small R/^ 



In order to get an analytical result we use the property 
of products of Besscl functions [13, HE] 

J,+„{A)JniB) = - fde {A Bey e-»^ 

(A2) 
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with X"^ = + — 2AB cos{9), which in our case reads 



Jl{Rk) 



d9Jo 



2i?/Cfcsin 



(A3) 



We introduce this in Eq. (jAl[) and interchange the inte- 
grations 



Jo 



(A4) 



Then we are able to calculate the inner integral 
Jo [4i?/esin(f)y] 



m 



2 / dy- 

'0 



(A5) 



= 2/n 



2i?/esin( - )y 



2R/^sm[-]y 



with the modified Bessel functions of zeroth order Iq and 
Kq . In order to simplify the integral we split it at a small 
cutoff e := 2arcsin((5), with S = {2R/^Ci)-^ small for 
sufficiently large R/S,. Here Ci is an unknown numerical 
constant. The cutoff e,S ^ 1 is chosen such, that one 



can use the expansion of Io{z)Ko{z) « 1/(22;) (for large 
argument z) at > e. The equation for the fluctuations 
simplifies to 



+ / dOe-'" h 
Jo 




-ie 



4i?/esin(f) 



(A6) 



2i?/^sin( - 



2i?/^sm \-]y 



One can show that the second integral vanishes for e <C 
1, in spite of the logarithmic divergence of Kf) at small 
arguments. Carrying out the first integral for e, (5 <C 1 
we finally obtain the particle number fluctuations 



{5N^) ^^\n [CR/^] , 
9 



(A7) 



with C :— 2Cie~^+'"^. By numerically calculating the 
integral we find C ~ 7.4. This corresponds to Ci ~ 13.7 
and thus 5 k, (27i?/^)^^. Hence e,(5 ^ 1 is satisfied as 
long as i?/<^ > 0.04. Fig. (O compares Eq. (|X7)) with 
the numerical solution of the integral in Eq. (|Al[) . 
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